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ABSTRACT 

The ubiquity of filamentary structure at various scales through out the Galaxy has triggered a 
renewed interest in their formation, evolution, and role in star formation. The largest filaments 
can reach up to Galactic scale as part of the spiral arm structure. However, such large scale 
filaments are hard to identify systematically due to limitations in identifying methodology 
(i.e., as extinction features). We present a new approach to directly search for the largest, 
coldest, and densest filaments in the Galaxy, making use of sensitive Herschel Hi-GAL data 
complemented by spectral line cubes. We present a sample of the 9 most Herschel 

filaments, including 6 identified from a pilot search field plus 3 from outside the field. These 
filaments measure 37-99 pc long and 0.6-3.0 pc wide with masses (0.5-8.3)x 10^ Mq, and 
beam-averaged (28", or 0.4-0.7 pc) peak H 2 column densities of (1.7-9.3)x 10^^ cm“^. The 
bulk of the filaments are relatively cold (17-21 K), while some local clumps have a dust 
temperature up to 25-47 K. All the filaments are located within <60 pc from the Galactic mid¬ 
plane. Comparing the filaments to a recent spiral arm model incorporating the latest parallax 
measurements, we find that 7/9 of them reside within arms, but most are close to arm edges. 
These filaments are comparable in length to the Galactic scale height and therefore are not 
simply part of a grander turbulent cascade. 
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1 INTRODUCTION 

The Milky Way is a barred spiral galaxy. Stars form in dense molec¬ 
ular clouds concentrated on the Galactic disk. Giant molecular 
clouds (GMCs), the birth places of high-mass stars, are often or¬ 
ganized in complex filamentary networks (e.g.. Bally et al. 1987; 
Busquet et al. 2013). Extremely long and filamentary GMCs can 
reach up to Galactic scale as part of spiral arm structures. Exam¬ 
ples are the two most prominent infrared dark clouds (IRDCs): the 
“Nessie” and “Snake” nebulae. Goodman et al. (2014) present an 
updated view of the physical Galactic plane and a careful analysis 
of CO data and find that the 80 pc IR-dark Nessie (Jackson et al. 
2010) is part of a much longer bone-like structure (up to 5 times 
the IR-dark extent) tracing the center of the Scutum-Centaurus 
spiral arm. On the northern part of the same arm lies the Snake, 
a >30 pc long sinuous IR-dark filament, although its molecu¬ 
lar extent has not been explored (Wang et al. 2014; Wang 2015). 
The Nessie and Snake both contain a chain of ~1 pc clumps 
likely fragmented from the pristine filament; some of those mas- 
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sive (10^ — 10^ M©) clumps are collapsing and forming high-mass 
star clusters (Jackson et al. 2010; Henning et al. 2010; Wang et al. 
2014; Wang 2015). Thus, Snake- and Nessie-like filaments are di¬ 
rectly linked to high-mass star/cluster formation, and bridge local 
to Galactic structures. These large scale filaments are one to two 
orders of magnitude larger than the filaments found in the Gould 
Belt clouds (Andre et al. 2014). 

If similar large scale filaments are common features of spiral 
arms, many should be observed towards other lines-of-sight. How¬ 
ever, the Snake and Nessie are the only two examples reported to 
be associated with the Milky Way’s spiral arm structure. Hardly 
any other similar filaments are known (cf. Ragan et al. 2014) de¬ 
spite the fact that the largest molecular cloud complexes are indeed 
good tracers of spiral structure (Dame et al. 1986). Moreover, pre¬ 
vious studies associate filaments with spiral arms in the longitude- 
velocity space alone. A recent spiral arm model has incorporated 
parallax distance of the arms (Reid et al. 2014), allowing a direct 
association with important distance information (§ 4.2). 

Very long and skinny IRDCs like the Snake and Nessie are 
rare features in IR surveys (Peretto & Euller 2009) because (a) they 
need a favorable location against background in order to be seen 
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Figure 1. Three steps in finding large scale filaments, (a) Composite image where 500/350/70 /rm images are color coded as red/green/blue. The red ellipse 
outlines the candidate filaments, (b) Column density (cm“^, upper) and dust temperature (K, lower) maps, (c) Customized PV diagram made following the 
curvature of the filaments. Clearly two filaments (G28, G29) are shown in the lower panel. The red line in the upper panel denotes the PV slice path, and the 
red cross marks the overlap of the two filaments in 2D. 


in extinction, and (b) they are usually made of individual compact 
extinction regions that are not always obviously connected. This 
background dependent method is therefore biased (as reflected by 
the term “infrared dark cloud” itself). See Wilcock et al. (2012) and 
Ragan et al. (2014) for more discussion on the limitations in the 
extinction method. 

A more robust way to identify filaments is to directly see them 
in emission. Since the bulk of the filaments are cold (< 20 K), they 
are most readily detected in far-IR. In this study, we make use of the 
Herschel Infrared Galactic Plane Survey (Hi-GAL, Molinari et al. 
2010) to conduct the first systematic search for Snake- and Nessie- 
like filaments in the Galaxy. We demonstrate the feasibility of our 
new identification procedure and present first results from the lon¬ 
gitude range of 15° < I < 56°. In forthcoming papers we will 
extend the study to the entire Galactic plane complemented with 
ongoing spectral line surveys and investigate structure and kine¬ 
matics of the identified filaments. 


2 DATA 

The Herschel key project Hi-GAL (Molinari et al. 2010), the first 
unbiased Galactic plane survey at far-IR wavelengths, provides a 
sensitive data set to directly search for large scale, cold and dense 
filaments. Hi-GAL covers the entire Galactic plane with nominal 
\b\ ^ 1° (following the Galactic warp). The data set is comprised of 
continuum images at 70, 160, 250, 350, and 500 /im obtained with 
PACS (Poglitsch et al. 2010) and SPIRE (Griffin et al. 2010) cam¬ 
eras onboard the Herschel Space Observatory (Pilbratt et al. 2010), 
with nominal beam sizes of 5.2", 12", 18", 25", and 37", respec¬ 
tively. The flux uncertainty in Hi-GAL data is typically 20%. 

Spectral line cubes are used to check the velocity coherence of 
the filaments and verify that a given 2D filament is indeed a single 
object, not multiple, unrelated objects along the line of sight. We 
retrieve the (1-0) data cubes from the Galactic Ring Survey 
(GRS, Jackson et al. 2006). The GRS covers the northern Galactic 
plane (18° ^ I ^ 55.7°, |5| ^ 1°; and 15° ^ < 18°,|5| ^ 

0.5°), with a resolution of 46" and an RMS noise of 0.4 K per 
0.2 km s“^ channel. 

Our sample also includes three filaments outside the GRS cov¬ 
erage: the Snake (Gll), Nessie (G339), and G64 (§4.1). For these 
filaments we use additional spectral line data. We mapped the en¬ 
tire extent of Gll in and (2-1) using the Atacama 

Pathfinder Experiment (APEX) 12m telescope. The observations 


were conducted on 2013 September 4, 7, 8, and 9, in four on-the- 
fly maps. The full width at half maximum (FWHM) beam is ~ 28" 
and the la rms is 0.2 K per 0.1 km s“^ channel. For G339 we use 
^^CO(l-O) from the Three-mm Ultimate Mopra Milky Way Sur¬ 
vey (ThrUMMS; Nguyen et al. 2015)^. For G64, we use the CfA 
CO (1-0) survey (Dame et al. 2001). Detailed analysis of the spec¬ 
tral line data will be presented in a forthcoming paper (Wang et al. 
in preparation). 


3 METHODS 

The structure of the interstellar medium is hierarchical and fila¬ 
mentary at multiple spatial scales. Depending on the decomposi¬ 
tion method, filaments can be very complex (e.g. Schisano et al. 
2014). Here, we focus on the largest, coldest, and densest filaments 
with simple morphology seen in Hi-GAL images. We identify cold 
filaments with three criteria: morphology, temperature, and veloc¬ 
ity coherence, as illustrated in Fig. 1. (a) Morphology. Using the 
Herschel images of the “Snake” and “Nessie” as a guide, we have 
searched for similar features by visual inspection. A filament is de¬ 
fined by a skinny long feature with aspect ratio ^10 and showing 
high contrast with respect to its surroundings. As we are looking for 
the early phases in evolution, we use the long wavelength images 
at 350 and 500 /xm primarily, (b) Temperature. We then construct 
temperature and column density maps from pixel-by-pixel SED fit¬ 
ting (Appendix B). Filaments that exhibit systematically lower tem¬ 
perature with respect to their surroundings are selected. The mor¬ 
phology is verified using column density maps, (c) Velocity coher¬ 
ence. Finally, all filaments are checked for velocity coherence in a 
customized position-velocity (PV) diagram/6>//6>wmg the curvature 
of the filament. Coherence means continuous, not broken, emission 
and velocity as a function of position along the filament. 

Compared to the traditional extinction approach, this direct 
procedure removes geometric ambiguity and bias. One important 
advantage is the temperature information, enabling us to select the 
cold (thus pristine) clouds, ideal for studying the onset of star for¬ 
mation. For each selected filament, we remove the Galactic scale 
background emission from the Herschel images using a Fourier 
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Transform based routine (Appendix A). We then fit the spectral en¬ 
ergy distribution (SED) to the multi-wavelength images to obtain 
temperature and column density maps (Appendix B). The codes 
are made publicly available (see Appendixes). 


4 RESULTS AND DISCUSSION 

4.1 A sample of the 9 most prominent filaments 

We have applied the identification procedure to the GRS field 
and identified 6 filaments. Outside the GRS field, we include also 
Snake, Nessie, and G64 which is identified by chance. They are 
of the same nature as others and we include them in an effort to 
compile the most prominent filaments known to date. In total, we 
present 9 filaments (Table 1, Fig. 2). Note that the selected sam¬ 
ple is not meant to be complete; instead, it should be regarded as 
the first descriptive sample of the most prominent, cold and dense, 
large scale filaments like Nessie and the Snake in the Galaxy. They 
illustrate the most important visual features and physical character¬ 
istics needed to classify GMCs as filaments. 

Table 1 lists physical parameters of the sample. The caption 
describes how the parameters are derived. Overall, the filaments 
measure 37-99 pc long and 0.6-3.0 pc FWHM wide at distances of 
2.8-5.4 kpc, with masses (0.5-8.3)x 10^ M©, and beam averaged 
peak H 2 column densities of (1.7-9.3)x 10^^ cm“^. The bulk of 
the filaments are cold (17-21 K), while some clumps have higher 
temperature of 25-47 K most plausibly due to local star formation 
activities. The filaments show remarkable velocity coherence, with 
small mean velocity gradients of 0.07-0.16 km s“^ pc“^ measured 
along the filaments. In certain clumps, however, the velocity gradi¬ 
ents are much larger. The linear mass densities along the filaments 
are 10^ — 10^ M© pc“ ^. The maximum linear mass density that can 
be supported by thermal pressure (Wang et al. 2011, 2014), given 
the lowest temperature within each filament, lies in the range of 
28-34 M© pc“^. If the filaments are made of bundles of purely 
thermally supported sub-filaments like the Taurus L1495/B213 fil¬ 
ament (Hacar et al. 2013), each of these filaments must be com¬ 
prised of 4-36 thermal sub-filaments. It is more plausible that they 
are made of turbulence supported sub-filaments, because the large 
velocity dispersion of 1.4-3.1 km s“^ suggests the presence of 
supersonic turbulence spread throughout these filaments, similar 
to high-mass star forming regions and IRDCs (Wang et al. 2009, 
2011,2012, 2014). 

We classify the morphology of these filaments into five cat¬ 
egories: L: linear straight or L-shape; C: bent C-shape; S: quasi- 
sinusoidal shape; X: crossing of multiple filaments; H: head-tail or 
hub-filament system. A “head” is a prominent/bright clump resides 
at the tip of a filament. Some filaments are characterized by more 
than one class (Table 1). For example, G49 is head-tail (H), and 
another filament crosses at the middle of the tail (X). These fila¬ 
ments represent the simplest morphological modes of filamentary 
GMCs. Different morphology may have resulted from different fil¬ 
ament formation processes (e.g., Myers 2009). Understanding the 
physics of these simple modes is the first step towards understand¬ 
ing more complex filamentary systems. While smaller filaments, 
from sub-pc to 10 pc scales, naturally occur in turbulent boxes (e.g. 
Moeckel & Burkert 2014), the filaments we observe are too large 
to be part of this process. 

100 pc appears to be the upper limit of these Herschel iden¬ 
tified filaments, whose minimum column density is larger than 
~10^^ cm“^. Some, if not all, of them are the dense part of larger. 


Large-scale filaments in Milky Way 3 



Figure 2. Gallery of the selected filaments: column density maps plotted in 
logarithmic scale in Galactic coordinates. See peak/mean densities for each 
filament in Table 1. The scale bars denote 10 pc. The red dotted line depicts 
Galactic plane at the source distance; when accompanied by a red arrow, it 
means that the plane is outside of the plotted area, parallel to the red dotted 
line, at the direction of the arrow. Table 1 lists distance to Galactic plane for 
each filament. 

less dense structures. For example, Nessie can be traced up to 430 
pc in CO (Goodman et al. 2014). Another example is G26, which is 
the dense and cold part of a much larger CO filament (Ragan et al. 
2014). The largest coherent CO filament reported so far is the 500 
pc “wisp” located at I ~ 50° (Li et al. 2013). Those “CO filaments” 
have minimum column densities ~10^° cm“^ (Ragan et al. 2014), 
an order of magnitude lower than our ''Herschel filaments”. 


4.2 Galactic distribution 

All the filaments are located within Galactic mid-plane, with height 
“above/below” the mid-plane |^|<60 pc. For comparison, the 
FWHM thickness of the molecular disk at the Galactocentric radii 
of these filaments (4-8 kpc) is 80-100 pc (the Galactic scale height 
is FWHM/V81n2 ~ 40 pc) (Nakanishi & Sofue 2006), and the in¬ 
trinsic spiral arm widths are 170-630 pc (Reid et al. 2014). These 
filaments are comparable to or larger than the Galactic scale height 
and therefore are not simply part of a grander turbulent cascade. 

The exact Galactic distribution of these filaments depends on 
the adopted spiral arm model. In Fig. 3 we compare the location 
of the filaments to the latest spiral arm model of Reid et al. (2014). 
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Table 1. Physical Parameters of the Selected Sample 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

ID 

Name 

"^LSR? '^max 

(7 V 

Distance 

Size 

Tdust 

mean/max 

Mass 

z 

Class 

Arm 



(km s“^) 

(km s“^) 

(kpc) 

(pcxpc) 

(K) 

(10^^ cm“^) 

(104 Mq) 

(pc) 



Gll 

CFGOll.11-0.12 

[27, 30.5, 33] 

1.9 

S 4Q+*^-^2 

'^•^^-0.36 

37x0.6 

17-25 

0.8/2.3 

1.1 

11 

S 


G24 

CFG024.00±0.48 

[93, 96.0, 99] 

3.1 


82x1.9 

17-31 

0.5/2.6 

8.1 

59 

C 


G26 

CFG026.38±0.79 

[44, 47.7, 50] 

1.4 

q 1 q-|-0.24 
'^•^'^-0.25 

61x1.1 

18-25 

0.3/2.1 

1.9 

62 

X 

Scu. 

G28 

CFG028.68-0.28 

[84, 88.2, 92] 

3.0 

4.89+°il 

60x1.6 

17-32 

0.8/7.6 

5.5 

-8 

s 

Scu. 

G29 

CFG029.18-0.34 

[89, 93.8, 98] 

3.5 


99x1.6 

20-28 

0.5/2.0 

5.5 

-15 

L 

Scu. 

G47 

CFG047.06±0.26 

[53, 57.5,61] 

2.3 


78x3.0 

17-29 

0.2/1.5 

2.0 

39 

C 

Sag. 

G49 

CFG049.21-0.34 

[66, 68.5, 72] 

1.9 

5 41 +0-31 
^•^-'--0.28 

85x1.6 

21-40 

0.9/7.5 

8.3 

-13 

H,X 

Sag. 

G64 

CFG064.27-0.42 

[18, 22.0, 25] 

2.0 

'^•^^-1.56 

51x2.7 

17-35 

0.1/1.7 

0.5 

-3 

L,H 

Loc. 

G339 

CFG338.47-0.43 

[-43,-37.5, -35] 

3.7 

Q oq-|-0.26 

Z.OO_o 28 

80x1.0 

17-47 

0.6/9.3 

2.9 

-2 

S,H 

Scu. 


Col. (1): ID number used throughout this paper. Col. (2): Approximate centroid Galactic coordinates in degrees with heading “CF” (cold filament). Col. (3-4): 
Velocities and velocity dispersion determined from ^^CO(l-O), with exception of Gil [^^CO(2-l), Wang et al. in prep] and G64 [CO (1-0), Dame et al. 
(2001)]. The Umin and Umax are determined from a PV plot customized following the curvature of the given filament (Fig. 1); ulsr and velocity dispersion 
(7v are determined from Gaussian fitting to a spectrum averaged over the entire filament. Col. (5): Kinematic distance calculated using the prescription of 
Reid et al. (2009) with the latest Galactic parameters from Reid et al. (2014). For sources with distance ambiguity, near distance is adopted. The error corre¬ 
sponds to a quoted uncertainty of 5 km s“^ in systemic velocity. For G49 we adopt the parallax distance of W51 Main/South (Sato et al. 2010), which is very 
close to the derived kinematic distance for the filament (5.5 kpc). Col. (6): Projected length (following the filament’s curvature) and FWHM width (Appendix 
C), deconvolved with the 28" beam size. Col. (7-9): Range of dust temperature, mean/peak beam averaged column density, and mass measured from a polygon 
encompassing the filament. Col. (10): Height from the coordinates given in col. (2) to the physical Galactic plane taking into account (a) zq = 25 pc and (b) the 
dynamic center of the Milky Way at Sagittarius A* is / = 359.944°, h = —0.046° (Reid & Brunthaler 2004; Goodman et al. 2014). Col. (11): Morphology 
class, see § 4.1 for details. Col. (12): Association with spiral arm according to the Reid et al. (2014) model (Fig. 3). 

Individual note: Gll — a well studied IRDC nick named “Snake” (Wang et al. 2014). G26 — dense part of a larger, less dense filament reported by 
Ragan et al. (2014). G28 — reported by Schneider et al. (2014) but without discussion in the context of Galactic filaments. Associated with IRDC G28.53—0.25 
(Swift 2009). G49 — associated with a well studied Hit region/GMC complex W51. The filament associated with W51 has been reported previously 
(Carpenter & Sanders 1998; Parsons et al. 2012; Ginsburg et al. 2015). Parameters for G49 are derived for the filament only, excluding the Hit region which is 
at a slightly different i^lsr- G339 — IRDC also known as “Nessie” (Jackson et al. 2010; Goodman et al. 2014). 



X (kpc) 


Figure 3. Distribution of the 18 Galactic scale filaments known to date 
plotted on the latest Milky Way spiral arm model (Reid et al. 2014). The 
white circular sector highlights our search field. Solid curves are the cen¬ 
ters of the best fitted arms and the corresponding dotted lines show their 
widths. Distance error bars for the ''Herschel filaments” from this work are 
from Table 1; for “CO filaments” from the literature (Beuther et al. 2011; 
Battersby & Bally 2012; Li et al. 2013; Tackenberg et al. 2013; Ragan et al. 
2014) a uniform error bar of ±0.5 kpc is plotted for reference. 


The model was derived from VLBA parallax measurements of over 
100 masers and represents the most precise model so far. We find 
that 7 out of the 9 (78%) Herschel filaments (red symbols in Fig. 3) 


are associated with spiral arms within their distance error, whereas 
6 of them reside near arm edges and only 1 (G49) close to an arm 
center. Here we have included G339 (Nessie) as “associated” if we 
extrapolate the Scutum arm model (which is limited by observed 
data points in Reidet al. 2014). On the other hand, Gll (Snake) 
and G24 are outside the Icr width of the arms. This is surprising 
because Nessie and the Snake are located right in the center of 
the Scutum arm when examined in the longitude-velocity plot of 
the Vallee (2008) and Dame & Thaddeus (2011) spiral arm models 
(Wang et al. 2014; Goodman et al. 2014). Expanding the statistics 
to all the 18 large scale filaments known to date (red and blue sym¬ 
bols in Fig. 3), 12 of them are associated with arms (extrapolation 
of the arm model is needed for two filaments), whereas only 2 are 
close to an arm center. 

If we compare to the Vallee (2008) spiral arm model in the 
traditional fashion, i.e., in / — u space (e.g. Ragan et al. 2014; 
Goodman et al. 2014), 5/9 of the Herschel filaments, or 8 of the 
18 Herschel and CO filaments, are associated with arms. The only 
filament in our sample that is outside arms in both models is G24 
(however it does not stand out in any of the physical parameters 
listed in Table 1, cf. Eden et al. 2013). The inconsistency is due to 
the difference in the two spiral arm models. Because the Reid et al. 
(2014) model has accurate distance, it should be regarded as the 
most up-to-date model. In summary, comparing to the Reid et al. 
model, 67% (12/18) of all known large scale filaments are located 
within spiral arms, whereas 11% (2/18) reside close to an arm cen¬ 
ter. A simple statistical test shows that our Herschel filaments lie 
preferentially on spiral arms, while the CO filaments identified by 
Ragan et al. (2014) is randomly distributed (Appendix D). 

How many large scale filaments are there in the Galaxy? Our 
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systematic search provides a rough estimate. In the GRS field (/ 
runs from 15° to 56°) we have identified 6 filaments. The largest 
distance is close to kpc. This search field is a circular sector 
that occupies 12.88 kpc^ of the Galactic plane area (Fig. 3). If large 
scale filaments are distributed uniformly in the disk at Galactocen- 
tric radii of i? G [3,8.5] kpc, we estimate about 90 Nessie-like 
features. Similarly, Ragan et al. (2014) identified 7 CO filaments in 
the GRS field where the most distant one is at ^4 kpc, implying 
about 240 similar CO filaments within the solar circle. 

In the coming years, we plan to expand the search to the entire 
longitude range using Hi-GAL and ongoing spectral line surveys. 
The ongoing BeSSeL project (Reid et al. 2014) and in the near fu¬ 
ture, GAIA, will continue to improve the measurement of spiral 
arms. Spatially resolved simulations have started to provide quan¬ 
titative constraints on the distribution of large scale filaments with 
respect to galaxy discs and spiral arms (Smith et al. 2014; Dobbs 
2015). Fully operational ALMA and NORM A will be capable of 
resolving similar filaments in the LMC, SMC, and nearby galax¬ 
ies (for reference, 1 arcsec at 10 Mpc corresponds to 48 pc) (e.g. 
Schinnerer et al. 2013; Pety et al. 2013). Synergizing this informa¬ 
tion will eventually reveal a full, hierarchical picture of Galactic 
star formation in the context of filamentary structure. 
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Figure Al. Example of FT based background removal, performed on a 
2° X 2° region centered on the filament of interest (red box). 


APPENDIX A: FT BASED BACKGROUND REMOVAL 

For proper analysis of the structure within the filaments, one has 
to remove emission from the background and foreground. In the 
literature, different methods have been used to define background: 
interpolation across masked-out regions (Peretto et al. 2010), Gaus¬ 
sian fitting along Galactic latitude (Battersby et al. 2011), or using 
a constant value from a reference region (Juvela et al. 2012). These 
methods involve assumptions about the structure of the background 
or do not allow for spatial variation in background. We use a new 
method based on Fourier Transform (FT). We transform the orig¬ 
inal image into Fourier domain and separate low and high spatial 
frequency components, then inverse FT the two components sep¬ 
arately back into the image domain. Low spatial frequency com¬ 
ponent corresponds to large scale background/foreground emis¬ 
sion, while high spatial frequency component corresponds to small 
scale structures. Fig. Al shows an example in a field approximately 
2° X 2° centered on the filament G28. One can see that the method 
separates the Galactic background very well. 

The separation between low and high spatial frequencies is 
determined from the power spectrum. In these images, most power 
comes from large scale emission, corresponding to low spatial fre¬ 
quencies (the central pixels in power spectrum, equivalent to “short 
spacing” in the (u, v) sampling of interferometric observations). 
Keeping this monotonic relation in mind, we start with an initial 
threshold of 90% of the maximum of the power spectrum. Pixels 
above this cut are masked (selected), setting other pixels into zero 
power, and are inverse Fourier transformed to produce the low spa¬ 
tial frequency map. Similarly, pixels below the cut are transformed 
to produce the high spatial frequency map. Images such as Fig. Al 
are made and inspected, then the cut is changed by a step of ±0.5%. 
This process is iterated until the background image becomes rep¬ 
resentative of the Galactic scale variation. Our experience shows 
that ~90% is a good value for most images. The structure in the 
resulting background emission has a scale of >0.5°. 

The idea of this method originates from a well-known “short 
spacing” problem in radio interferometric imaging, where the cen¬ 
tral hole in the (u, v) space caused by a lack of the shortest 
baselines leads to large scale emission being filtered out. Similar 
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FT-based filtering is also widely used under the name of “unsharp 
masking” in graphic processing. The method is particularly use¬ 
ful in determining Galactic scale background in Herschel images, 
but can also be applied to images of any angular size. Compared 
to other methods, this approach does not involve any assumption 
about the underlying structure and picks up spatial variation in 
background. On top of that, the method is simple and easy to im¬ 
plement. We make the routine publicly available ^ . 

After removal of background, the derived column density is 
up to 20-30% lower and the dust temperature is 0.2-2 K lower for 
the bulk, cold region of the filament, than derived from the original 
images. 


APPENDIX B: SED FITTING 

To convert the Herschel continuum images to physical parameters, 
we have developed a procedure to fit a modified black body func¬ 
tion to the multi-wavelength images on a pixel-by-pixel basis ^ . In 
this study, images at 70, 160, 250 and 350 /xm are used while the 
500 /xm image is excluded because of its low resolution. All im¬ 
ages are first convolved to a circular Gaussian beam with FWHM 
= 2S" using the kernels provided by Aniano et al. (2011), and then 
re-gridded to the same pixelization (8'' pixel size). The convolu¬ 
tion takes into account the small asymmetry in the Herschel beams. 
Pixels below the median value in the 350 /xm images are masked 
out, and this mask is applied to other wavelengths. This mask out¬ 
lines the approximate boundary of the filament. For each unmasked 
pixel, intensity at various wavelengths are modeled as 

= B.(l-e-"-), 

where Planck function Bu is modified by optical depth 
(Kauffmann et al. 2008) 

Tv — 

which is caused only by dust. Here, gas to dust ratio is as¬ 
sumed to be i^gd = 100. Dust opacity per unit dust mass is 
(Ossenkopf & Henning 1994) 

Dust emissivity index is fixed to /3 = 1.75 in the fitting. The free 
parameters are the dust temperature and column density expressed 
in number of molecular hydrogen per unit area. The fitting provides 
Tdust and A^H 2 maps along with error maps for further analysis. 
Typical error in Tdust is <2 K and <20% in iVH 2 • 


APPENDIX C: CROSS-FILAMENT PROFILES 

Density and temperature profiles are extracted from 1-3 cuts per¬ 
pendicular to a given filament. The cuts are chosen to be wide 
enough to represent the general properties of the filament, while 
narrow enough to avoid curvature. In cases where one cut cannot 
represent the width of the entire filament, multiple cuts are placed 
in different representative segments. Averaged profiles from the 
cuts are used for analysis. Fig. Cl shows a representative exam¬ 
ple. As one can see, the density and temperature profiles are anti¬ 
correlated, i.e., denser regions are in general colder. The FWHM 

^ https://github.com/esoPanda/FTbg 

^ The SED fitter has been implemented as a standalone python package 

(http : / /hi-gal - sed- fitter . readthedocs . org). 



Figure Cl. Temperature and density profiles of G64. The cut is perpendic¬ 
ular to the filament; the profiles are averaged within the thickness of the cut 
(17.30 which is chosen to cover the cold bulk (the central three clumps) of 
the filament. 

widths of the filaments are measured by fitting a Gaussian to the 
density profile after removal of the “emission free” baseline (cut 
<2 and >9 pc in Fig. Cl), and are listed in Table 1. Uncertainty 
of the width is 10-30% depending on how many cuts are averaged. 
In a forthcoming paper, the profiles will be analyzed with detailed 
physical models to investigate the structure and stability. 


APPENDIX D: STATISTICAL TEST 

In the same search field, we identified 6 ''Herschel filaments” with 
5 (83%) associated with arms, while Ragan et al. (2014) identified 
7 “CO filaments” with 3 (43%) associated with arms. We regard 
a filament to be associated with arm if it lies within the width of 
the arms within its distance error (Fig. 3). Given the large distance 
error bars, however, the random chance to be associated with arm 
is relatively high towards some lines of sight. 

We perform a simple statistical experiment to test the random 
probability. We draw a series of line of sight rays every 2° of longi¬ 
tude within the GRS field, and for each ray we calculate the prob¬ 
ability of a random filament along this line of sight to fall within a 
spiral arm. The probability depends on the distance error. For error 
of 0, ±0.25 (appropriate for our Herschel filaments within the GRS 
field), and ±0.5 kpc (appropriate for the Ragan et al. 2014 CO fil¬ 
aments), the probability averaged for all the lines of sight is 25%, 
37%, and 48%, respectively. 

Comparing these random probabilities to the observed frac¬ 
tions, the test demonstrates that the Herschel filaments do prefer¬ 
entially lie on spiral arms (83% vs 37% random), and the on-arm 
fraction of the Ragan et al. (2014) filaments is most likely random 
(43% vs 48% random). This is consistent with the fact that the Her¬ 
schel filaments have higher average column densities. 
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